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This study determined tine individual maps of all fingers in Brodmann area 3b of the human 
primary somatosensory cortex in a single fMRI session by tactile stimulation at 19 sites 
across all phalanges and digit bases of the 5 right-hand digits. To quantify basic features 
of the digit maps within and across subjects, we applied standard descriptive measures, 
but also implemented a novel quantitative analysis. This so-called Direction/Order (DiOr) 
method tested whether subjects exhibited an ordering of peak fMRI representations 
along their individual direction of alignment through the set of analyzed phalanges and 
whether these individual directions were similar across subjects. Across-digit analysis 
demonstrated that for each set of homologous phalanges, the D5-to-D1 representations 
were successively represented along a common direction of alignment. Hence, the well- 
known mediolateral D5-to-D1 somatotopy was not only confirmed for the distal phalanges 
(pi ), but could also be shown for the medial (p2) and proximal phalanges (p3). In contrast, 
the peak activation for the digit bases (p4) only partly elicited that digit succession. 
Complementary, intra-digit analysis revealed a divergent picture of map topography for the 
different digits. Within D5 (and in a trend: D4), an ordered p1-to-p3 succession was found 
across subjects, pointing to a consistent intra-digit somatotopy for D5, with p3 generally 
found medial-posterior to pi. In contrast, for D1 , D2, and D3, most subjects did not present 
with ordered p1-to-p3 maps nor were directions of alignment similarly oriented between 
subjects. These digits therefore exhibited highly diverse representation patterns across 
subjects. 
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INTRODUCTION 

Primary sensory cortices are characterized by topographic maps 
which represent the peripheral receptor sheet in a highly ordered 
fashion (Kaas, 1997). In the auditory cortex, frequency is encoded 
along a single tonotopic axis, whereas the retinotopic map in the 
primary visual cortex incorporates two perpendicular visual-field 
axes to encode eccentricity and polar angle. The somatotopy of 
the human primary somatosensory cortex (SI), first described 
by Penfield's somatosensory homunculus (Penfield and Boldrey, 
1937; Penfield and Rasmussen, 1950), comprises the tactile body- 
surface representation along the postcentral gyrus, representing 
the digits D5 (little finger) to Dl (thumb) in medial- to-lateral 
succession. Electrophysiological mapping of the hand area in 
anesthetized primates implies a second axis perpendicular to 
the first, representing the phalanges within a digit, with the 
fingertip or distal phalanx (pi) being represented toward the 
Brodmann area (BA) 3a/3b border and the proximal parts of 
the digit toward the BA 3b/ 1 border (for new- world owl mon- 
keys, see Merzenich et al., 1978; for old-world macaques, see Kaas 
et al., 1979 and Nelson et al., 1980). Electrophysiological map- 
ping data from awake monkeys is sparse, but in awake macaques, 
the intra-digit representations of the phalanges were found to be 
less ordered (Iwamura et al, 1983). Neurons with receptive fields 
for the proximal phalanges were still found close to the BA 3b/ 1 



border, whereas representations of distal phalanges were not lim- 
ited toward BA 3a but present across the entire respective digit's 
BA 3b strip. 

In humans, several studies confirmed the medial-to-lateral 
succession of the digit representations for D5 through Dl in 
SI, with functional magnetic resonance imaging (fMRI) during 
tactile or electrical stimulation delivered to the distal phalanx 
(e.g.. Nelson and Chen, 2008; Schweizer et al, 2008; Sanchez- 
Panchuelo et al., 2010). So far, similar across-digit somatotopies 
of the second or third phalanges have not been investigated in 
humans by imaging. Attempts to obtain intra-digit maps, by stim- 
ulating the phalanges within a digit, have so far mostly been 
limited to single digits and yielded inconsistent results both with 
fMRI and MEG (magnetoencephalography). An overview of all 
previous studies on intra-digit mapping in individual humans 
is given in Table 1. In an fMRI study, Blankenburg et al. (2003) 
stimulated the three phalanges (pi to p3) and the digit base (p4) 
of D3. Combining the data of eight subjects into a group anal- 
ysis, they found a succession of pl-to-p4 activations within BA 
3b, pointing toward the BA 3b/ 1 border. This pattern is similar to 
the intra-digit somatotopy obtained in anesthetized owl monkeys 
(Merzenich et al, 1978). Focusing on the left hand, Sanchez- 
Panchuelo et al. (2012) applied the traveling-wave stimulation 
and analysis approach (originally developed for mapping the 
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Table 1 | Overview of fMRI and MEG studies of intra-digit topology on individual human subjects. 



Publication 


Technique 


Subjects* 


Hand 


Digits 


Phalanges 


Cortical region 


Reported intra-digit pattern 


Blankenburg 
et a!., 2003 


fMRI 


6 


Right 

(dominant) 


D3 


pi to p3** 


BA 3b and BA 1 


Ordered pi to p3 representations ortlnogonal to 
central sulcus (inferior to superior, statistically 
significant) 


Overduin and 
Servos, 2004 


fMRI 


6 


Right 

(dominant) 


D1 
D2 
D4 


pi to p3 


BA 3a, BA 3b, 
and BA 1 


Only individual results, across-subject 
somatotopy not evaluated 



Schweisfurth fMRI 10 Right D2 pi and p4 BA 3b No consistent topology between pi and p4 

etal., 2011 (dominant) representations across subjects 

D5 Representation of p4 medial to pi (statistically 

significant) 



Sanchez- 
Panchuelo 
etal., 2012 



fMRI 



Left 



D2 pi to p4 BA 3b, 1, 2 Ordered pi to p3 representations orthogonal to 

central sulcus, with reversals between BAs, in 
4 out of 6 subjects 



Sanchez- 
Panchuelo 
etal., 2014 



fMRI 



Left D2 to D4 pi to p4 

simultaneously 



BA 3b, 1, 2 Ordered pi to p3 representations orthogonal to 
central sulcus, with reversals between BAs, in 
3 out of 4 subjects 



Current study fMRI 11 Right D1 pi to p3 BA 3b No consistent topology between pi , p2, and p3 

(dominant) representations across subjects 

14 D2 No consistent topology between pi , p2, and p3 

representations across subjects 
12 D3 No consistent topology between pi , p2, and p3 

representations across subjects 
11 D4 Ordered p1-p3 representations along course of 

CS (anterior to posterior, statistical trend) 
11 D5 Ordered p1-p3 representations along course of 

CS (lateral-anterior to medial-posterior, 

statistically significant) 



Hashimoto 
et al., 1999a 



Hashimoto 
et al., 1999b 



MEG 



MEG 



14 



Right 



Right 



D2 pi to p3 BA 3b 

(plus three 
points 
at palm) 

D2 pi to p3 BA 3b 

(plus two 
points 
at palm) 



No consistent topology between equivalent 
current dipole (ECD) pi , p2, and p3 
representations across subjects 



No consistent topology between equivalent 
current dipole (ECD) pi , p2, and p3 
representations across subjects 



Hlushchuk MEG 11 Right D2 pi and p3 BA 3b Representation of p3 superior to pi (statistically 

et al., 2004 (dominant) significant) 

Tanosaki and MEG 17 Right D3 pi and p3 BA 3b Representation of p3 lateral to pi (statistically 

Hashimoto, significant) 

2004 

'Only subjects who showed significant activation in respective area. 

''Stimulation involved pi to p4 and palm. Because only one subject showed significant activation for all five locations, only pi to p3 were taken into account for 
this table. 



visual cortex, Sereno et al, 1994) to the phalanges and base of D2. 
They showed individual activation patterns in which they visu- 
ally identified a mirror- reversal pattern for pi to p3 and p3 to pi 
across SI in 4 out of 6 subjects, suggesting an ordering similar as 
in anesthetized macaques (Nelson et al, 1980). The group further 
showed similar phase maps in 3 out of 4 subjects jointly for 
the digits D2 to D4, by simultaneously stimulating the phalanges 



pi to p3 (sequentially) of these three digits (Sanchez-Panchuelo 
et al., 2014). In 2004, Overduin and Servos explored individ- 
ual intra-digit phase bands for Dl, D2, and D4 in BA 3b, but 
stated that they did not investigate across-subjects somatotopy 
due to difficulties in map interpretation caused by limitations in 
cortical surface sampling. Also analyzing individual maps elicited 
by tactile stimulation of the fingertip (pi) and the base (p4) of 
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the little finger (D5) and the index finger (D2), our group could 
show that in BA 3b the base presentation was consistently located 
medial to the fingertip representation for D5, whereas for D2 no 
consistent across-subject pattern was observed between the rep- 
resentations of the fingertip and the base (Schweisfurth et al., 
2011). 

Using MEG and stimulating pi to p4 of the index finger, 
Hashimoto et al. ( 1999a,b) repeatedly observed non-ordered pha- 
lanx representations, whereas Hlushchuk et al. (2004) found p3 
to be represented superior to pi, in qualitative agreement with 
the monkey data. Exploring the phalanx representations of the 
middle finger, Tanosaki and Hashimoto (2004) localized the p3 
representation lateral to that of pi, in contrast to the Blankenburg 
et al. (2003) where p3 was found superior to pi. 

To gain further insight into the individual intra-digit maps of 
specific digits, we used fMRI to determine the BA 3b representa- 
tions for all phalanges of all five digits of the right, dominant hand 
at an isotropic spatial resolution of (1.5 mm) ^. Optimization of 
the stimulation paradigm enabled us to measure all phalanges and 
the digit base of each finger in one fMRI measurement and the 
data from all 19 stimulation sites within one fMRI session. Based 
on these data the general layout as well as the individual variations 
of the BA 3b maps of the entire finger area could be investigated 
across a larger cohort of healthy subjects. 

The large number of stimulation sites also called for a more 
quantitative approach in describing the spatial pattern of phalanx 
representations than just visual inspection or manual delineation. 
With the Euclidean distance measure we specifically analyzed the 
spatial locations of the peak activations in each subject along two 
axes: First, along the 5 digits (across-digit map), where the dis- 
tances between the representations of the homologous phalanges 
(distal phalanges, second phalanges, third phalanges, digit base) 
were explored separately across digits, and secondly, within indi- 
vidual digits (intra-digit map), considering the distances between 
the different phalanges (pi to p4). 

To quantify the main directions along which the representa- 
tions were oriented and to compare these individual BOLD acti- 
vation patterns across subjects, we used a novel analysis approach 
termed Direction/Order (DiOr) method. Principle component 
analyses were applied to the peak-activation locations for either 
the across-digit maps or the intra-digit maps. Hereby, we could 
identify the main direction along which the homologous pha- 
langes across digits or the different phalanges within a digit were 
distributed in each subject and compare these individual direc- 
tions statistically across subjects. In addition to this main axis 
along which the phalanx presentations were distributed, we also 
assessed whether activations associated with neighboring skin 
positions were ordered along the subject's main direction. If 
the main directions were similar across subjects and individual 
patterns were generally ordered along the subjects' main direc- 
tions, we assumed consistent across-subject somatotopy for the 
respective phalanx or digit. 

MATERIALS AND METHODS 
MRI 

Eighteen healthy subjects (6 women, range 21-30 years, 
mean 27 ± 2.5 years) participated in two MRI sessions. 



The study was approved by the ethics committee of the 
Georg-Elias-Miiller-Institute of Psychology, Gottingen 
University, and informed written consent was obtained before 
each session. All subjects were right-handed according to the 
Edinburgh Inventory (laterality index 0.9 ±0.1, Oldfield, 1971). 
Magnetic resonance imaging (MRI) was performed at 3 T (TIM 
Trio, Siemens Healthcare, Erlangen, Germany) using a 32- 
channel head coil. The first session incorporated a whole-brain 
structural MRI measurement, a partial-volume localizer fMRI 
measurement for delineating the presumed digit-representation 
area along the central sulcus, whole-brain echo-planar imaging 
(EPI), and MR angiography. The second session consisted of 
5 fMRI measurements, each mapping the phalanges within a 
single digit. 

Structural MRI was based on sagittal Tl-weighted 3D 
MPRAGE (magnetization-prepared rapid gradient-echo) acqui- 
sition [repetition time (TR) = 2530 ms, echo time (TE) = 

3.4 ms, flip angle = 7°, acquisition matrix = 256 x 256, 160-192 
partitions, resolution = 1 x 1 x Imm^, total acquisition time 
(TA) = 10:49 min]. The images were used to identify the 
motor hand knob at the pre-central gyrus (Yousry et al, 1997) 
and to perform surface reconstructions of the cortical white- 
matter/gray-matter boundary. 

Functional images were recorded using a gradient-echo EPI 
sequence at 1.5 x 1.5 x 1.5 mm^ resolution (TR = 2000 ms, 
TE = 36ms, flip angle = 70°, acquisition matrbc = 128 x 128, 
field of view = 192 x 192 mm^, partial Fourier factor = 6/8). 
Nineteen double-oblique (transverse-to-sagittal and transverse- 
to-coronal) sections were positioned perpendicular to the central 
sulcus, cutting the motor hand knob in mediolateral direc- 
tion and covering its whole depth. The AutoAlign Scout feature 
(Siemens Healthcare, Erlangen, Germany) ensured an identical 
slice position for all consecutive fMRI runs for the individual 
subject. 

A single gradient-echo EPI volume was performed with iden- 
tical orientation and position as the functional measurements, 
but covering the entire brain with 81 sections (1.5 xl.5x 

1.5 mm^ resolution, TR = 8600 ms, TE = 36 ms, flip angle = 
70°, acquisition matrix = 128 x 128, field of view = 192 x 
192 mm^, partial Fourier factor = 6/8), to facilitate registra- 
tion of the partial-volume images to the whole-volume structural 
images. MR angiography was based on a Tl-weighted 3D FLASH 
sequence (TR = 22ms, TR = 4.43ms, flip angle = 18°, resolu- 
tion = 0.3 xO.3 X 0.5 mm^, 57 sections from two overlapping 
slabs) with identical volume coverage and slice orientation as the 
functional images. 

For three of the subjects, the second MRI session with all- 
phalanx mapping was repeated to assess the retest-reliability of 
the resulting maps by calculating the Euclidean distance between 
the peak-vertex coordinates of the two mapping sessions. 

TACTILE STIMULATION AND fMRI MAPPING 

Tactile stimulation was applied by a piezo-electric unit 
(QuaeroSys, St. Johann, Germany) controlling independent 
stimulation modules, each with an 8-dot Braille display (2x4 
matrix, covering an area of 2.5 x 7.5 mm^) at its distal top 
face (Figure lA) (Schweizer et al, 2008; Schweisfurth et al.. 



Frontiers in Human Neuroscience 



www.frontiersin.org 



September 2014 | Volume 8 | Article 658 | 3 



Schweisfurth et al. 



Comprehensive digit map of liuman BA 3b 




■ I iniiii I iiiiii 

0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 

Time / images 

FIGURE 1 I Tactile stimulation. (A) Piezo-electric device used for 
vibrotactile stimulation of individual phalanges. (B) Stimulated digit areas for 
mapping (pi , red; p2, blue; p3, green; and p4, yellow). (C) One digit was 
stimulated per run, the example refers to the positioning of the piezo-electric 
devices (black frames in B) to the ring finger. (D) Pseudo-randomized 
stimulation design (color code as above, baseline = gray). 



2011). Stimulation frequency was set to 32 Hz. In each cycle 
two out of the 8 pins were randomly chosen and elevated 
(square wave, stimulation duration = 10.4 ms, inter-stimulus 
interval = 20.8 ms) resulting in a fast varying stimulation pattern 
across the entire Braille display. The pins were set to maximum 
drive-out, to obtain a salient tactile stimulation which was clearly 
perceivable. Throughout the measurements subjects were told to 
keep their right hand relaxed and pronated. 

For the digit localizer measurement, five stimulation modules 
were individually positioned under the distal phalanx (pi) of each 
digit of the right hand, with the Braille display being oriented cen- 
trally below the phalanx, parallel to the digit axis, with the bulk 
of the module pointing away from the digit. During the fMRI 
measurement sequential stimulation was applied, separately stim- 
ulating each fingertip from the thumb (Dl) to the little finger 
(D5), with a rest period embedded after each fourth digit stim- 
ulation (Dl, D2, D3, D4, Rest, D5, Dl, D2, D3, Rest, D4, D5, 
Dl, D2. . . ). The duration of the stimulation at each site and for 
each rest period was 12 s (6 images). During the entire measure- 
ment each of the 5 digits was stimulated 8 times (total = 96 s, 
48 images). Total measurement time including 10 distributed rest 
periods (total = 120 s, 60 images) and baseline (20 s, 10 images) 
added up to 10:20 min (310 images). 

For intra-digit phalanx mapping tactile stimulation was 
applied to the first, second, and third phalanx (pi, p2, p3) and 
to the base (p4) of each digit (Figure IB). The term "base" refers 
to the volar pads described in primate studies (Merzenich et al., 
1978) and denotes the palmar skin position over the caput of the 
metacarpal bone proximal to the respective digit (Blankenburg 
et al., 2003). Each of the five fMRI measurements in the second 
MRI session explored the topography of the phalanges within a 



single digit. Four stimulation modules (only 3 for the Dl due 
to the missing p3) were individually positioned centrally below 
the phalanx, but orthogonally to the digit axis (Figure IC). For 
stimulation of the phalanges of D2 and D3 the bulk of the mod- 
ules pointed to the ulnar, for D4 and D5 to the radial, and for 
Dl to the ulnar (digit base and p2) and distal (pi) side. Each 
functional measurement (Figure ID) consisted of 6 cycles of a 
stimulation block (average: 84 s, 42 images) alternating with a rest 
block (average: 12 s, 6 images). Within this stimulation block, the 
distal phalanx (pi) was stimulated once (12 s on average: jitter 
between 8, 12, or 16 s) and each of the other phalanges (p2, p3) 
and the base (p4) twice (2 x 12 s on average: jitter between 8, 
12, or 16 s) within the block to account for the lower sensitiv- 
ity of p2, p3, and p4 (Johansson and Vallbo, 1979; Schweisfurth 
et al, 2011). The succession of the four stimulation sites within 
the stimulation block was pseudo-randomized such that stimula- 
tion of one phalanx was always followed by stimulation of another 
phalanx or a rest block. During fMRI of one digit, the total stim- 
ulation time was 72 s (36 images) for the first phalanx (pi) and 
144 s (72 images) for the phalanges p2, p3, and p4. Total scan 
time for the functional measurement, including baseline (20 s, 
10 images), was 9:56 min or 298 images (7:32 min, 226 images 
forDl). 

To assure attention on the stimulation, each functional mea- 
surement was associated with a task in which subjects were 
instructed to covertly count short randomly distributed inter- 
rupts embedded in the tactile stimulation (duration 156 ms, 
repetition every 0.5-3 s). After completion of the measurement, 
subjects reported the overall number of interrupts they counted 
and received verbal feedback from the experimenter about the 
actual number of interrupts during the measurement. 

PRE-PROCESSING AND CO-REGISTRATION 

The first part of data analysis was performed with Brain Voyager 
QX 2.3 (Goebel et al., 2006, Brain Innovation, Maastricht, The 
Netherlands). Structural Tl-weighted images were corrected for 
intensity inhomogeneities, the axial slice orientation was brought 
parallel to the plane through the anterior and posterior commis- 
sures (AC-PC plane), and a cortical mesh at the white-matter 
gray-matter border was reconstructed. Functional runs were 
motion-corrected in k-space (Siemens Healthcare, Erlangen, 
Germany) and image space (Brain Voyager) as well as temporally 
high-pass filtered as implemented by default in Brain Voyager. In 
the same step, each functional run for each subject was registered 
to the first fMRI measurement (functional localizer) of the first 
session (6 degrees of freedom, using trilinear/sinc interpolation) 
to ensure an accurate alignment between all functional measure- 
ments within a subject. This first functional measurement of each 
subject was registered (6 degrees of freedom) to the individ- 
ual whole-brain Tl-weighted structural image; additional manual 
optimization ensured its best possible registration specifically in 
the region of the central sulcus. Using trilinear interpolation, 
the resulting registration matrix was then applied to each func- 
tional run of the subject to achieve identical projections onto the 
Tl-weighted 3D data in AC-PC-oriented native space. During 
registration the fMRI data were also interpolated from 1.5 mm 
isotropic resolution to 1 mm isotropic resolution. No stereotactic 
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transformation was applied in any of the performed registration 
steps. 

The MR-angiography volume was evaluated within FSL 
(FMRIB Software Library, Smith et al, 2004), where it was reg- 
istered to the Tl -weighted image using FLIRT (FMRIB's Linear 
Image Registration Tool, Jenkinson and Smith, 2001). More pre- 
cisely, the registration of the MR-angiography volume to the 
Tl -weighted volume was achieved through concatenation of 
intermediate registration steps from the MR-angiography vol- 
ume to the first volume of the functional localizer measurement 
to the whole-brain EPI volume to the Tl -weighted volume. 
This resulted in transformation of the MR-angiography volume 
onto an interpolated 0.5 mm isotropic Tl-weighted volume in 
Brain Voyager AC-PC-oriented native space. 

The next step involved the exclusion of voxels around promi- 
nent blood vessels. This method is commonly applied in high- 
resolution fMRI studies of visual processing in humans (Cheng 
et al, 2001; Yacoub et al., 2001; Shmuel et al, 2007) to reduce the 
BOLD signal from these vessels, which limit the spatial specificity 
of the BOLD response (see Polimeni et al., 2010 for the visual sys- 
tem). For the somatosensory system this was based on a two-step 
approach. Using the functional localizer, fMRI activation clusters 
in which the peak vertices (of an activation cluster) of three or 
more digits were in very close proximity (not further than 3 mm 
apart in any direction) were identified, based on the assumption 
that such an extensive spatial coincidence of peak activation from 
the distal phalanx of three or more digits is not in line with the 
expected separate digit representations in BA 3b and could there- 
fore represent a BOLD signal change due to a large blood vessel 
(Schweizer et al., 2010). In those cases the respective location in 
the high-resolution angiograms was inspected for the presence 
of a vessel. If a vessel could be observed at the respective posi- 
tion, the area activated (at q(FDR) < 0.05) upon stimulation of 
the two assumingly most distally represented fingers (as expected 
from mediolateral fingertip somatotopy, e.g., D4 and D5 for posi- 
tions close to the expected Dl representation or Dl and D2 for 
positions close to the expected representations for D4 and D5) 
was marked and respective vertices were excluded from further 
analyses. 

DATA ANALYSIS 

Stimulus-related activations were determined by a GLM analy- 
sis. For the functional localizer measurement, the main effect for 
the predictor of each of the 5 distal phalanges was calculated. 
For the intra-digit measurements, a GLM analysis for each digit 
was conducted calculating the main effects for the 4 predictors 
pi to p4 (only pi, p2, and p4 for the thumb). All predictors were 
convolved with a two-gamma function (Friston et al., 1998) and 
voxels passing a positive f-value threshold corresponding to a false 
discovery rate of q(FDR) = 0.05 were considered significant. The 
resulting maps were projected onto the individual white-matter 
gray-matter surface reconstruction. The value at each vertex was 
calculated by taking into account the values of the statistical map 
voxels up to 3 mm orthogonal to the mesh, as it is standard 
in Brain Voyager 2.3. Further analysis of the resulting maps was 
restricted to BA 3b, here defined as being located at the poste- 
rior wall of the central sulcus. Because no in-vivo marker for the 



identification of BA 3b borders in individual subjects is avail- 
able, somatosensory mapping studies use a schematic approach 
to define the putative BA 3b/ 1 border (suggested by Moore et al., 
2000, used e.g., by Blankenburg et al, 2003; Nelson and Chen, 
2008, and Stringer et al., 2011). Here, the junction of the supe- 
rior portion of the posterior bank of the central sulcus with the 
beginning of the curve of the crest of the post-central gyrus was 
individually defined as the putative border between BA 3b and 
BA 1 (Moore et al, 2000). The medial-to-lateral extent of the 
putative BA 3b within this border was based on the extend of 
the individual pi activations of the 5 digits measured in the first 
functional measurement, plus an addition of approximately 1 cm 
to either side of the most medial and lateral activation. After BA 
3b definition, the mesh-spanning vertices of this area and their 
thresholded t-values within this BA 3b digit area were exported 
to MATLAB. 

For further analysis, the coordinates of the mesh vertex with 
the highest f-value were determined for each phalanx of each fin- 
ger. This peak- value approach is generally applied in somatosen- 
sory mapping studies (Kurth et al., 2000; Blankenburg et al., 2003; 
Nelson and Chen, 2008; Schweizer et al., 2008; Stringer et al., 
2011). It does not take into consideration the description and 
comparison of the BOLD-activation overlap, which is beyond the 
scope of the present study. 

To provide a distance measure for comparison of the present 
data with other studies, mean Euclidean distances were calculated 
for each phalanx within a set of homologous phalanges, between 
the phalanx' peak vertex and the respective thumb's phalanx peak 
vertex: 

/ {Xdigit.phal — XD\,phal )^ + {y digit, phal. — yDl.phal.)^ 
y + iZdigit,phal. - ZDl,phal.)^ 

Analogously, for each phalanx within a digit, the mean intra- 
digit Euclidean distance between the phalanx peak vertex and the 
respective digit's first-phalanx peak vertex was calculated by 

/ (^digit ,phal. ~ Xdigit,pl)^ + {ydigit,phal. ~ y digit, plf' 
y + iZdigit,phal. ~ Zdigit,pl)^ 

To describe features of the map outline both within and across 
subjects, across-digit and intra-digit analyses were conducted, 
using the newly-developed two-step DiOr method as described 
below in detail for across-digit analysis. 

In a first step of across-digit analysis, respective activation 
peaks were analyzed to determine the direction of largest varia- 
tion for each set of homologous phalanges (pi across 5 digits, p2 
across 5 digits, p3 across 4 digits, p4 across 5 digits). A principle 
component analysis (PCA) was applied to the peak-vertex coordi- 
nates of the phalanx representations in order to extract their main 
direction of variation, i.e., the direction in AC-PC-oriented native 
space along which they were mostly distributed. For each sub- 
ject and set of homologous phalanges the PC vector N = (Nx, Ny, 
Nz) was obtained (being the eigenvector belonging to the largest 
eigenvalue of the covariance matrix, normalized to a length of 1). 
It described the best linear fit and explained the largest amount 
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of variance across the peak positions within a set of homolo- 
gous phalanges (Figure 2). To ensure that the N vector pointed 
approximately into the direction of D1-D5, its sign was chosen 
such that the N coordinate with the largest absolute value had 
the same sign as the respective coordinate of the difference vector 
D5-D1 for across-digit analysis (i.e., if Nx>Ny and Nx>Nz, then 
if D5x-Dlx>0, also Nx > 0). To investigate the possibility of a 
common direction pattern across subjects, for each set of homol- 
ogous phalanges (and digit bases) the respective N coordinates 
along each of the three axes were compared across subjects to 
explore a possible across-digit direction consistency. To that end, 
first Fisher's transformation (take the atanh of the value, Fisher, 
1915) was applied to the individual N vector coordinates and then 
one-sample f-tests (two-tailed, atp < 0.017 corrected for testing 
along the three axes) against zero were calculated. 

In addition to the direction along which the representations 
were aligned, we explored in a second step whether a Dl-to-D5 
succession could be detected along that direction. For each sub- 
ject and set of homologous phalanges, it was assessed whether 
a Dl-to-D5 ordering was observed across the 5 digits along the 
individual direction of largest variation. For that purpose, the 
respective peak vertices were transferred into an orthonormal 
coordinate system which included N as first base vector. After 
this transformation, the first coordinate by definition described 
the position of each phalanx representation along the line given 
by N. Using this information, it could be assessed whether the 
peak vertices were ordered (from Dl to D5) according to the 
definition that neighboring phalanges yield similar or sequen- 
tially increasing coordinates along the PC vector and that not all 



peak vertices are identical. Across subjects, binomial tests (p < 
0.05) then examined whether more subjects than expected by 
chance showed a Dl-to-D5 succession for homologous phalanges 
(across-digit ordering). The DiOr approach used the probabili- 
ties to find such an ordered pattern in individual subjects (i.e., 
probabilities 0.056 for pi, p2, and p4 as well as 0.189 for p3, by 
assuming equal likelihoods for each possible pattern). 

For the intra-digit analysis, we similarly assessed the individual 
directions along which the representations of a digit's phalanges 
were aligned and possibly ordered. To that end, we performed 
analogous steps as in the across-digit analysis but including the 
peak-vertex coordinates of the three phalanges within a digit 
(two phalanges for the thumb). In each intra-digit analysis, only 
those subjects were included that showed significant activation 
for all phalanges of the respective digit. For each subject and 
digit, the normalized PC vector N was calculated, defined such 
that it approximately pointed in the direction of pl-p3 (from 
pi to p2 for the thumb). At group level, the N coordinates of a 
digit were compared to investigate a possible intra-digit direction 
consistency along any of the axes. For D2 to D5, a second analy- 
sis step assessed whether more subjects than expected by chance 
presented with ordered phalanges (from pi to p3, intra-digit 
ordering) along the PC vector (i.e., probability 0.5 for D2 to D5). 
For Dl, no such analysis was carried out, as there by definition 
always is a pl-to-p2 ordering. 

RESULTS 

In the present study, individual activation maps in BA 3b, cover- 
ing each of the phalanges and digit bases of all right-hand digits, 
were obtained and analyzed. Common features such as a similar 
main direction of succession across subjects and ordering along 
this direction were described on two levels: for homologous pha- 
langes (all pi, all p2, all p3, or all p4) across digits and for different 
phalanges within single digits. 

To control for attention on the tactile stimulation, subjects had 
to covertly count short randomly-distributed interrupts of the 
tactile stimulation. Because there was on average only a 4% dif- 
ference between the actual and reported numbers for each finger 
across subjects, it can be assumed that subjects were able to focus 
their attention toward the stimulation. 

Activation patterns acquired during the functional localizer 
measurements, stimulating the most distal phalanges (pi) of 
each of the five digits, generally showed the expected medial-to- 
lateral succession of digit representations from D5 to D 1 along the 
posterior wall of the central sulcus. 

After the first session, 3 of the 18 subjects had to be excluded, 
either due to no detectable activation in the assumed BA 3b, or 
dominance of a large vessel parallel to the entire SI digit area, or 
a rare variation of a divided central sulcus at the level of the hand 
knob (Alkadhi and KoUias, 2004). 

Out of the remaining 15 subjects, in a second session complete 
tactile stimulation of all 4 within-digit stimulation sites could be 
accomplished in all subjects for D5, in 14 subjects for D4, D3, and 
D2, and in 12 subjects for Dl (all 3 sites). Missing measurements 
resulted from difficulties in positioning the piezo-electric mod- 
ules or from fatigue of the subjects due to the lengthy session. 
After statistical analysis, projection to the surface reconstruction. 
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FIGURE 2 I Direction vector N (for pi of subject S2). (A) View from 
above onto tlie surface reconstruction of an individual left hemisphere, with 
the motor hand l<nob [Q.) in the foreground and the somatosensory BA 3b 
hand area in the center (A-P anterior-posterior; M-L, medial-lateral; l-S, 
inferior-superior). The area used for peak-vertex analysis is marked (white 
border), peak vertices are marked by white crosses. From left to right, digit 
activations refer to Dl, magenta; D2, yellow; D3, green; D4, blue; and D5, 
red. (B) Three-dimensional representation of the 5 peak vertices (1 mm 
isotropic coordinates in AC-PC-oriented native space) and their best-fitting 
line (dashed) as direction vector N. This subject shows somatotopy as the 
projections of the first phalanges onto N are sequential from Dl to D5. 
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and elimination of assumed large vessel contribution, significant 
BA 3b activations of all phalanges and the base within a digit were 
present in 11 subjects for D5, D4, and Dl, in 12 subjects for D3 
and in 14 subjects for D2. The complete map, i.e. detectable acti- 
vation from all phalanges and all bases of all digits (total of 19 
stimulation sites) was available in 7 subjects. 

Averaging across the 19 stimulation sites in these seven sub- 
jects, a significant activation consisted of 62.95 ± 86.51 voxels 
[mean ± standard deviation (SD)], and the average peak-voxel's 
t-value was 5.50 ± 1.59 (mean ± SD). Numbers of voxels and 
f-values per phalanx and digit are presented in Table 2. 

Visual inspection on the surface reconstruction confirmed 
the known succession of digit activations for the first phalanx 



(pi), progressing from medial-posterior-superior (D5) to 
lateral-anterior-inferior (Dl) (Figure 3, first column). The data 
from the additional stimulation sites revealed a similar pattern 
for the second and third phalanx, as both sets of homologous 
phalanges (all p2 and all p3) exhibited the analog succession 
from medial-posterior-superior (D5) to lateral-anterior-inferior 
(Dl) (Figure 3, second and third column). This succession is also 
evident in the peak- voxel distribution, as illustrated in Figure 4A, 
where the peak vertices of all stimulation sites are plotted in 3D 
coordinate space. The figure also shows that the peak phalanx 
activations of the same digit tend to group together, resulting in 
an overall D5-to-Dl succession of the phalanx representations 
across the five digits. The finger-base (p4) representations, in 



Table 2 | Volume and t-values of BOLD activations across and within digits. 



Volume: mean ± standard error of 
the mean (SEM) (median) [mm^] 



t-value: mean ± SEM (median) [mm^] 



Average across digit 


Average across phalanx 




Average across digit 


Average across phalanx 


Dl 48 ± 2 (25) 


pi 


82 ± 1 (59) 


Dl 


5.49 ± 0.05 (5.36) 


pi 


6.25 ± 0.03 (5.99) 


D2 58 ± 1 (32) 


p2 


69 ± 2 (47) 


D2 


5.26 ±0.03 (5.06) 


p2 


5.57 ± 0.03 (5.32) 


D3 62 ± 2 (40) 


p3 


47± 1 (24) 


D3 


5.67 ± 0.03 (5.30) 


p3 


5.27 ± 0.02 (5.06) 


D4 92 ± 3 (40) 


p4 


52 ± 2 (25) 


D4 


5.68 ±0.04 (5.30) 


p4 


4.76 ± 0.03 (4.75) 


D5 53 ± 1 (42) 






D5 


5.45 ± 0.03 (5.17) 







7776 first and third column present the average volume and t-vaiues across all phalanges for each digit. The second and fourth column depict the average volume 
and t-values across all homologous phalanges. 




FIGURE 3 I Across-dlgit representations of phalanges and bases 
(subjects S3, S5, and S7). View from above onto the inflated surface 
reconstruction of individual left hemispheres with the motor hand knob (Q) in 
the foreground and the somatosensory BA 3b hand area in the center 
(convex areas, light gray; concave areas, dark gray; putative vessels, white; 



A-R anterior-posterior; M-L, medial-lateral). The area used for peak-vertex 
analysis is marked (white border). Each column shows the activations for 
homologous phalanges pi , p2, and p3 as well as digit bases p4. Peak voxels 
are illustrated in the respective color. Dl, magenta; D2, yellow; D3, green; 
D4, blue, and D5, red. 
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FIGURE 4 I Organization of digit area in BA 3b. For eacli of the 7 subjects 
that presented with significant activation for all phalanges and bases, all 
peak vertices were normalized by the respective subject's mean value 
across phalanges (but not bases; A-R anterior-posterior; M-L, medial-lateral, 
l-S, inferior-superior; Dl, magenta; D2, yellow; D3, green; D4, blue; D5, red; 
all units in 1 mm). (A) Normalized peak-vertex coordinates for pi (=1), p2 
(=2), and p3 (=3). (B) Normalized peak-vertex coordinates for p4 (=4). 



contrast, showed a much more variable pattern across digits 
and subjects (Figure 4B), ranging from the vicinity of (or even 
overlap with) the respective phalanx activations (Figure 3, 
subject S3) to positions lateral to Dl (rarely) or more commonly 
medial to D5 (Figure 3, subject S7). 

This visually-observed succession for each of the three pha- 
langes and variability of the finger-base representations were 
reflected in the Euclidean distances between the peak vertices 
of homologous-phalanx representations (Table 3). The distance 
between Dl and D5 was 16.7 ± 4.2 mm for pi, 14.7 ± 7.0 mm for 
p2, 18.0 ± 3.2 mm for p3, and 11.2 ± 8.4 mm for p4. These rather 
similar D1-D5 distances for pi, p2, and p3 resulted from the 
strong overlap of the BOLD activations and the relatively closely 
located peak vertices of pi, p2, and p3 (and sometimes p4) within 
the digits. Compared to the D1-D5 distances for pi, p2, and p3, 
the value for the digit base (p4) tended to be smaller and more 
variable, pointing to the visually-observed variability of the p4 
peaks for the digits D3, D2, and Dl. 

The within-digit Euclidean distance between the proximal (p3) 
and the distal phalanx (pi) (Table 3) exhibited comparable values 
for the digits Dl through D4: 4.8 ± 2.6 mm for Dl, 6.2 ± 2.1 mm 
for D2, 6.6 ± 5.2 mm for D3, 5.0 ±3.1 mm for D4, and a reduced 
value of 3.4 ± 2.3 mm for D5. Here, standard deviations across 
measurements were comparable, pointing to a smaller pl-p3 dis- 
tance for D5. The Euclidean pl-p4 distance increased across the 
digits from D5 to Dl (Table 3), which has to be at least partly 
ascribed to the much higher spatial variability of the p4 repre- 
sentations of Dl, D2, and D3 being distributed across the whole 
digit area. 

Visual inspection of the within-digit maps (Figure 5) indi- 
cated that subjects generally showed a similar pl-to-p3 (or even 
pl-to-p4) succession within the expected digit area for D5 and to 
a lower degree for D4. AH other within-digit maps visually pre- 
sented with a high level of individuality, both at the volume and 
at the surface level. 

The retest reliability across sessions (evaluated in 3 subjects) 
revealed a mean Euclidean difference between the peak vertices 
of the two measurements of 2.9 ± 3.6 mm (averaged across the 
3 subjects and stimulation sites). 

In the following two subsections, we report the results of the 
statistical evaluation of the across- and intra-digit maps. 



Table 3 | Euclidean distance across digits and within digits. 



Distance to D1/mm 


D2 


D3 


D4 


D5 


pi (to pi of Dl) 
p2(to p2 of Dl) 
p3 (to p2 of Dl) 
p4 (to p4 of Dl) 


6.8 ±2.1 
4.6 ± 3.0 
8.1 ± 2.3 

6.9 ± 76 


11.4 ±3.2 
9.2 ± 5.1 
12.8 ±3.8 
12.6 ±9.7 


15.1 ± 5.7 
13.3 ± 70 
15.9 ±4.1 
8.2 ± 70 


16.7 ±4.2 
14.7 ± 70 
18.0 ±3.2 
11.2 ± 8.4 


Distance to p1/mm 


P2 


p3 


p4 




Dl (to pi of Dl) 
D2 (to pi of D2) 
D3 (to pi of D3) 
D4 (to pi of D4) 
D5 (to pi of D5) 


4.8 ± 2.6 
4.0 ± 2.8 

2.9 ± 2.7 
2.9 ±2.9 
2.3 ± 2.4 


6.2 ± 2.1 
6.6 ± 5.2 
5.0 ±3.1 
3.4 ± 2.3 


14.4 ± 8.7 
9.2 ± 4.2 
8.8 ± 5.6 
6.8 ±3.0 
5.4 ± 5.0 





For each phalanx within a set of homologous phalanges (upper part of table), 
the across-digit Euclidean distance (mean ± standard deviation) between the 
phalanx peak voxel and the respective thumb's phalanx peak voxel are stated 
(including the 7 subjects showing activation for each finger and phalanx). For 
each phalanx within a digit (lower part of table), the intra-digit Euclidean distance 
(mean ± standard deviation) between the phalanx peak voxel and the respective 
digit's first-phalanx peak voxel are stated (same subjects as above). 

STATISTICAL EVALUATION OF ACROSS-DIGIT MAPS WITH DiOr 

To statistically compare across-digit representation patterns 
across subjects, we determined the main spatial directions along 
which the homologous phalanges and the digit bases were aligned 
across the five digits (Table 4). For the distal, second, and third 
phalanges (pi, p2, and p3) individual direction vectors were 
similarly oriented across subjects, as averaging these individual, 
normalized vectors across subjects resulted in a mean direction 
vector of substantial and comparable length for the three pha- 
langes (see red vector in Figure 6A). Statistically, across-subjects 
t-tests (against zero) revealed significant across-digit direction 
consistency along all three axes, not only for pi but also for p2 and 
p3 (for p3 along the AP axis only in a trend; Table 4). The average 
direction vectors for the three phalanges pointed from lateral- 
anterior-inferior to medial-posterior-superior. For the digit bases 
(p4), such an across-digit direction consistency could not be 
observed, as the direction vectors of different subjects pointed to 
different directions, resulting in a small average direction vector 
across subjects (Figure 6A, Table 4). 

To determine a possible ordering between homologous pha- 
langes or digit bases along the direction vector, it was tested 
whether the projections of the respective peak vertices onto their 
best-fitting line (i.e., the direction vector N) were ordered from 
Dl to D5. For pi, p2, and p3, most subjects displayed this 
expected Dl-to-D5 ordering which could statistically be con- 
firmed not only for the first but also for the second and third pha- 
langes (Figure 6B). In the 6 out of 25 individual cases in which the 
expected across-digit ordering was not observed, the swap consis- 
tently occurred between neighboring digits: three times between 
D4/D5 (see e.g.. Figure 3, pi in Subject S5), twice between D3/D4 
(see e.g.. Figure 3, p3 in Subject S5), and once between D1/D2. 
For the digit base, on the contrary, no across-digit ordering 
could be observed, as only 2 out of 9 subjects showed Dl-to-D5 
ordering. 
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FIGURE 5 I Intra-digit representations of phalanges and bases 
(subjects S3, S5, and S7). View from above onto the inflated surface 
reconstruction of the left hemisphere with the motor hand knob (fi) in 
the foreground and the somatosensory BA 3b hand area in the center 
(convex areas, light gray; concave areas, dark gray; putative vessels. 



white; A-R anterior-posterior; M-L, medial-lateral). The area used for 
peak-vertex analysis is marked (white border). Each column shows the 
activations for phalanges pi (= red), p2 (= blue), and p3 (= green) as 
well as base p4 (= yellow) of a single digit (Dl to D5), peak voxels 
are illustrated in the respective color. 



Putting together the results from the direction-consistency 
analysis and the ordering analysis, for pi, p2, and p3 we generally 
found the Dl-to-D5 representations to be ordered along a simi- 
lar axis across subjects, pointing to an across-digit somatotopy for 
each of these phalanges. 

STATISTICAL EVALUATION OF INTRA-DIGIT MAPS WITH DiOr 

To statistically compare intra-digit representation patterns across 
subjects, for each subject and digit the main intra-digit direction 
was determined (Figure 7A), i.e., the principal direction along 
which the phalanges within a single digit were aligned in a sub- 
ject. Only for the digits D4 and D5, similar directions were found 
across subjects (Table 5, Figure 7A). D5 exhibited a significant 
intra-digit direction consistency from lateral-anterior to medial- 
posterior. For D4, an anterior-to-posterior intra-digit direction 
consistency was observed. The direction vectors of Dl, D2, and 
D3 did not show any consistency across subjects. 

To determine the intra-digit ordering for each digit (except for 
Dl), it was tested whether the projections of the peak vertices 
onto their direction vector were ordered from pi to p3. Only for 
D5 (but for no other digits) more subjects (9 out of 11) than 
expected by chance exhibited an ordered phalanx representation 
(Figure 7B). For D4, significance was missed by only a single 
subject, while no intra-digit ordering could be observed for D3 
and D2. 

Putting together the results from the direction-consistency 
analysis and the ordering analysis, for D5 we generally found the 
pl-to-p3 representations to be ordered along a similar axis across 
subjects, pointing to an intra-digit somatotopy for the little fin- 
ger. A similar trend was observed for the ring finger, while Dl, 
D2, and D3 presented with highly individual intra-digit maps. 



As additional test, pi to p3 of D3 and D5 were measured in six 
different subjects with a three-fold higher stimulation time per 
phalanx (unpublished data). For D5, the results showed a simi- 
lar N-vector as in the present study yielding pl-to-p3 ordering in 
5 out of 6 subjects. For D3, despite the increased amount of mea- 
surement time, again in line with the current work only 3 subjects 
presented with ordered pl-to-p3 representations. 

DISCUSSION 

BOLD MRI activation peaks, elicited by tactile stimulation of the 
phalanges and digit bases of all digits, were analyzed to quantita- 
tively describe the somatotopic maps in BA 3b within and across 
individual subjects. Due to the complexity of the data, they were 
not only described by standard measures as Euclidean distances, 
but also evaluated by the suggested DiOr approach, looking at the 
main direction of alignment and a possible ordering along that 
direction. Across-digit analysis considering homologous pha- 
langes across the 5 digits showed the Dl-to-D5 succession of 
the distal, medial, and proximal phalanges along the central sul- 
cus. The digit bases presented with a partially divergent cortical 
pattern. Intra-digit analysis revealed increasing Euclidean dis- 
tances along the phalanges, but only D5 and D4 presented with a 
consistent direction of alignment and only D5 with an ordering of 
pi, p2, p3 along this direction. Hence, only for D5 (and in a trend: 
D4) intra-digit maps resembled each other between subjects. 

COMPREHENSIVE FUNCTIONAL MAPS OF ALL PHALANGES AND DIGIT 
BASES 

The challenge to map the entire digit area of the right hand was 
met by combining an optimized stimulus paradigm including an 
attentional task with high-resolution fMRL The attempt allowed 
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Table 4 | Across-digit direction consistency. 

Phalanx p1 p2 p3 p4 

M-L A-P l-S M-L A-P l-S M-L A-P IS M-L A-P l-S 



Mean N -0.39 0.47 0.68 -0.28 0.52 0.70 -0.34 0.43 0.65 0.11 -0.06 0.01 

SD ±0.23 ±0.34 ±0.14 ±0.24 ±0.28 ±0.22 ±0.27 ±0.36 ±0.35 ±0.51 ±0.54 ±0.74 

p-value 0.001* 0.003* <0.001* 0.013* 0.002* 0.012* 0.012* 0.026 0.009* 0.520 0.770 0.466 

Mean N vectors (calculated across D l-to-D5 peak vertices), respective coordinates and standard deviations (SD) as well as p-values for phalanges pi to p3 and digit 
bases p4 along the medial-lateral (M-L), anterior-posterior (A-P), and Infehor-superior (l-S) ACPC axis, respectively. 
*p < 0.0167 (I.e., p < 0.05 corrected for the three comparisons along the three axes). 
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FIGURE 6 I Across-digit direction consistency and ordering. 

(A) Across-digit direction vectors for pi to p4 (black arrows, individual 
direction vectors; red arrows, mean values averaged across subjects; A-R 
anterior-posterior; M-L, medial-lateral; l-S, inferior-superior). Across-digit 
direction consistency is indicated by an asterisk and respective p-values. 

(B) Across-digit ordering as given by the number of subjects showing (black) 
or not showing (gray) ordered activations from D1 to D5 along the individual 
N vectors. The red lines illustrate the minimal number of ordered subjects 
needed to observe across-digit ordering. Phalanges showing ordering are 
indicated by an asterisk and respective p-values. 



for measuring all phalanges and digit bases in one session. Under 
these conditions, significant activations could be obtained for 
most stimulation sites in 15 subjects, out of which 7 subjects 
showed activation for all 19 stimulation sites. The sizes of the 
elicited activations were in the expected range. The similar sta- 
tistical f-values of the activations at the different stimulation 
sites approved our approach of stimulating the medial and distal 
phalanges as well as the digit base twice as often as the distal pha- 
langes, thereby compensating for the weaker fMRI signal due to 



the smaller cortical representations of the former. The exemplary 
re-test measurements demonstrated comparable maps across ses- 
sions, with a distance of less than 3 mm (<2 voxels) between 
respective peak vertices of the two sessions. 

ACROSS-DIGIT FINDINGS 

The quantitative activation maps of the homologous phalanges 
and digit bases across digits in BA 3b showed the succession of the 
phalanges from D5 to D 1 along the general course of the central 
sulcus (from medial to lateral, superior to inferior, and posterior 
to anterior). 

The increase of the average Euclidean distances between the 
peak representations of the distal phalanges of D2, D3, D4, and 
D5 and the one of Dl is generally interpreted as a reliable indi- 
cator of the expected succession of the individual digit represen- 
tations from D5 to Dl along the central sulcus (e.g.. Nelson and 
Chen, 2008) and is in complete concordance with previous fMRI 
studies describing the layout of the human digit area in BA 3b 
(Maldjian et al, 1999; Kurth et al, 2000; Nelson and Chen, 2008; 
Schweizer et al., 2008; Sanchez-Panchuelo et al., 2010). However, 
this interpretation can be problematic, as the Dl-to-D5 repre- 
sentations could still be distributed in an unordered way even if 
there is a succession in the distances. To conclude from increas- 
ing Euclidean distances that the representations are successive in 
space, all D2-to-D5 representations must in addition be found in 
the same direction of the Dl representation. This assumption is 
not required for the DiOr approach applied here, the results of 
which are discussed in the next paragraph. 

For the distal phalanx, the individual direction vectors (rep- 
resenting the main direction of alignment of the five activation 
peaks in each subject) all pointed into a similar direction. Further 
on, most of the subjects showed an ordered pi -representation 
succession from D5, D4, D3, D2 to Dl along this direction. These 
results demonstrate that the proposed DiOr analysis provides a 
valid and sensitive strategy to capture and statistically compare 
the layout of pi maps from individual subjects as well as to 
extract general features as the well-known medial-to-lateral pi 
somatotopy. 

The activations elicited by stimulation of the medial and 
proximal phalanges revealed an analogous pattern of increasing 
distances from Dl toward D2, D3, D4, and D5, hinting to an 
equivalent representation succession for the medial as well as for 
the proximal phalanges from D5 to Dl. This was further corrob- 
orated by the DiOr analysis, where both the medial as well as the 
proximal phalanges exhibited very similar across-digit directions 
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across subjects, again pointing along the central sulcus. Also, 
the representations were generally ordered along these directions, 
showing the corresponding succession from D5 to Dl for the 
medial phalanges and from D5 to D2 for the proximal phalanges. 
Consequently, not only the distal phalanges, so far used as rep- 
resentatives of the entire digits, showed this succession, but also 
the lower phalanges of the digits. Along with visual inspection, 
this finding indicates that the representations of all phalanges 
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FIGURE 7 I Intra-digit direction consistency and ordering. (A) Intra-digit 
direction vectors for Dl to D5 (black arrows, individual direction vectors; red 
arrows, mean values averaged across subjects; A-R anterior-posterior; M-L, 
medial-lateral; l-S, inferior-superior). Intra-digit direction consistency is 
indicated by an asterisk and respective p-values. (B) Intra-digit ordering as 
given by the number of subjects showing (black) or not showing (gray) 
ordered activations from pi to p3 along the individual N vectors. The red 
lines illustrate the minimal number of ordered subjects needed to observe 
intra-digit ordering. Digits showing ordering are indicated by an asterisk and 
respective p-values. 



of a single digit in BA 3b are located in relatively close vicinity, 
thereby forming a distinct area for each digit with relatively little 
overlap with the activations of other digits' phalanges. As these 
observations are analogous to the digit organization seen in elec- 
trophysiological mapping studies of squirrel, owl, and macaque 
monkeys (Merzenich et al, 1978; Sur et al, 1982; Iwamura et al., 
1983), this succession can be described as the axis of digit repre- 
sentation along the course of the central sulcus, in analogy to the 
eccentricity axis along the calcarine sulcus of the primary visual 
cortex (Sereno et al., 1994). We hence propose a well-ordered suc- 
cession of the Dl-to-D5 representations as single digit entities 
along the central sulcus, where each single digit representation 
area is composed of the representations of the three phalanges. 
Looking at individual subjects, Sanchez-Panchuelo et al. (2014) 
found such ordered digit entities for D2 to D4 by simultaneous 
stimulation of pi, p2, and p3. Our results are in line with and 
extend such data, as we could demonstrate a succession separately 
for each phalanx. 

For the digit bases (p4), we could not detect a consistent pat- 
tern of cortical representations across subjects. The Euclidean 
distance between the Dl representation and the D2, D3, D4, and 
D5 representations of the bases did not steadily increase such as 
in the phalanges. Digit-base peak representations were in some 
cases located close to the activations of the respective digit's pha- 
langes, but in other cases medial to the digit area of D5 or lateral 
to that of Dl. This variability was reflected both in the absence 
of a Dl-to-D5 ordering of the digit bases in most of the subjects 
and in the large between-subjects variation of the within-base 
direction vectors. Our observations resembled a mixture of sev- 
eral previous findings concerning the location of the digit-base 
representation. Penfield and Rasmussen (1950) distinguished a 
separate hand area, located medial to D5, from the representa- 
tional areas for each finger. In owl-monkey maps, the digit-base 
representations were described posterior to the pl-to-p3 repre- 
sentations (Merzenich et al., 1978). In macaques, the digit-base 
representations were found at both ends of the digit represen- 
tations: medial-posterior to the little finger for D5 and D4 and 
lateral-posterior to the thumb for D3 and D2 (Kaas et al, 1979; 
Nelson et al, 1980; Iwamura et al, 1983). 

INTRA-DIGIT FINDINGS 

The intra-digit analysis showed a much lower degree of consis- 
tency across subjects than the across-digit analysis. For each digit, 
the average Euclidean distance to the pi representation increased 



Table 5 | Intra-digit direction consistency. 
Digit D1 



M-L 



A-P 



l-S 



M-L 



D2 



A-P 



l-S 



M-L 



D3 
A-P 



l-S 



M-L 



D4 
A-P 



l-S 



M-L 



D5 
A-P 



l-S 



Mean N -0.27 0.16 -0.25 -0.12 0.06 -0.17 -0.10 0.21 -0.24 -0.30 0.37 -0.22 -0.43 0.46 0.23 
SD ±0.59 ±0.40 ±0.64 ±0.56 ±0.55 ±0.64 ±0.50 ±0.58 ±0.62 ±0.46 ±0.33 ±0.69 ±0.39 ±0.34 ±0.58 

p-value 0.137 0.238 0.280 0.526 0.940 0.360 0.822 0.449 0.303 0.134 0.004* 0.192 0.009* 0.003* 0.162 

Mean N vectors (calculated across the p7 to p3 peak vertices), respective coordinates and standard deviations (SD), and p-values for digits Dl to D5 along the 
medial-lateral (M-L), anterior-posterior (A-P), and Inferior-superior (IS) ACPC axis, respectively. 
*p = 0.0167 (I.e., p < 0.05 corrected for the three comparisons along the three axes). 
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from the medial to the proximal phalanx to the digit base. This 
might generally be interpreted as hint for a succession of the 
phalanx representations across the single-digit area. However, as 
discussed above, this only holds true if all lower phalanges are 
represented toward the same direction from the respective pi. As 
discussed in the next paragraph, this was only the case for some 
digits. 

Analyzing the individual directions along which the intra-digit 
representations were aligned and potentially ordered revealed a 
more differentiated view of the layout of the maps within the sin- 
gle digits. Across subjects, only the digits D5 and D4 presented 
with a similar axis along which phalanx activations were dis- 
tributed and only D5 exhibited an ordering from pi to p3 (for 
D4, significance was missed by one subject). 

In Table 1, all previous studies on intra-digit mapping in 
individual humans have been summarized. In the following para- 
graphs, the current findings are discussed in the light of this 
literature. 

For D5, our current observations confirm and extend the 
results of our previous study (Schweisfurth et al, 2011), demon- 
strating a consistent arrangement between pi and p4 for D5 
but not for D2 across subjects. The common direction along 
which the D5 representations were ordered pointed from lateral 
(pi) to medial (p4) in the previous and lateral- anterior (pi) to 
medial-posterior (p3) in the present study. Both results indicate 
an intra-digit succession along the general course of the cen- 
tral sulcus, not perpendicular to it as in anesthetized owl and 
macaque monkeys (Merzenich et al., 1978; Nelson et al, 1980). 
One reason for this finding may be the relatively small extent of 
the D5 digit area. This small D5 area is suggested by the observed 
Euclidean pl-p3 distances, which showed a tendency to decrease 
from Dl toward D5, a trend in accordance with the maps of awake 
macaques (Iwamura et al., 1983). The claim is also corroborated 
by inspection of other mapping studies, where the volume of the 
pi representation in BA 3b seems to decrease from Dl to D5 
and to get closer to the one of BA 1 (Kurth et al, 2000; Nelson 
and Chen, 2008). The smaller anatomical D5-area extent along 
with the small distance between the BA 3b D5 representation and 
the putative cytoarchitectonical BA 3b/ 1 border might explain the 
shift of the p2 and p3 representations into the medial direction. 

For D3, D2, and Dl, the directions along which the pha- 
langes were arranged differed substantially between subjects. In 
addition, no intra-digit ordering between phalanges along the 
individual direction vectors was found for D3 or D2. For D2, 
these results are in agreement with our previous study, where 
no consistent arrangement was observed between tip and base 
(Schweisfurth et al., 2011). However, this is in contrast to stud- 
ies in anesthetized monkeys (Merzenich et al., 1978; Nelson et al., 
1980) and two single-digit flVIRI studies (Blankenburg et al., 2003; 
Sanchez-Panchuelo et al., 2012). Blankenburg and colleagues 
reported a rostral(pl)-to-caudal(p4) succession of phalanx rep- 
resentations of D3 across BA 3b despite coarser spatial resolution 
and as a result of a group analysis of 8 subjects. Factors con- 
tributing to the divergence between their result and the present 
work may be our increased spatial resolution, our higher num- 
ber of subjects, and our analysis of individual maps. In the study 
of Sanchez-Panchuelo et al. (2012), the authors identified map 



reversals of the pl-to-p3 representations of D2 in SI, although 
no objective measure was given to justify the identification of 
the mirror-reversal positions. The across-subjects variation in the 
direction of the mirror reversal of the maps as well as the pres- 
ence of ordered maps in only 4 out of 6 subjects seem to be 
comparable to the variations in direction and ordering of our 
intra-digit representations of D2 (ordered maps observed in 8 
out of 14 subjects). The divergent results for anesthetized mon- 
keys also have to be put into perspective by the awake-macaque 
study (Iwamura et al, 1983) reporting that the narrow zone of 
proximal-phalanx (p2/p3) representations located toward the BA 
3b/ 1 border additionally contained receptive fields of distal pha- 
langes. This kind of intermittent spread of pi representations into 
the proximal-phalanx representation area could explain the diver- 
gent ordering and directions of phalanx representations within 
D3 and D2 among our subjects. 

Finally, the question arises whether methodological issues 
could be responsible for not observing consistent within-digit 
somatotopies across subjects in some digits. The presence of a 
consistent within-digit somatotopy in D5, exhibiting the smallest 
digit representational area of all digits, is a clear indicator that the 
data and the applied DiOr analysis were adequate to detect soma- 
totopic ordering across subjects. One limitation to be considered 
is the linearity of our DiOr approach, an anatomical simplifica- 
tion which does not take into account the complex anatomy of 
individual central sulci. We therefore carried out careful visual 
inspections of the Dl, D2, and D3 maps which confirmed the 
results of our objective analysis procedure and assured that no 
existing somatotopy was missed due to the linear simplification. 
The approach of measuring the complete map of all phalanges in 
one session (in order to omit the influence of across-session vari- 
ations) came with the price of a decreased number of averages per 
stimulation site, which in turn results in a slightly lower across- 
session reliability compared to our previous study (Schweisfurth 
et al, 2011). However, as results of additional measurements for 
D3 and D5 (unpublished data, see Results Section) with three- 
fold higher stimulation time are in line with the present results, 
the lower number of averages should not have been responsible 
for not observing somatotopy in some digits. 

The reported inter-individual map variations within some but 
not all digits correspond to the general context of individual 
variability of representational maps in the somatosensory cor- 
tex. Individual differences in digit maps could be explained by 
the individual usage of the fingers, shaping the cortical repre- 
sentational maps in a usage-dependent manner, as suggested and 
shown in monkeys (Merzenich et al, 1987; Jenkins et al., 1990; 
Recanzone et al., 1992) as well as in humans (Elbert et al., 1995; 
Sterr et al., 1998; Braun et al., 2000). The findings reported here 
might have resulted from a rather uniform use of the mainly "sup- 
porting" digits D4 and D5, contrasted by the individually varying 
use of the digits Dl to D3 of the dominant hand, which might 
have led to more individualized cortical intra-digit maps for these 
digits. 

CONCLUSIONS 

The present fMRI study describes the comprehensive map of 
all phalanges and digit bases across all digits of the dominant 
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hand. An in-depth exploration of individual digit somatotopy was 
performed both across and within all five digits in a larger num- 
ber of subjects. The novel quantitative approach DiOr described 
the individual directions along which the BA 3b BOLD-activation 
peaks aligned as well as the ordering along these directions. 

Separate analysis of the representations of the medial and 
proximal phalanges across the five digits revealed an ordered 
succession from medial (little finger) to lateral (thumb) along 
the central sulcus, extending the well-known Dl-to-D5 succes- 
sion of the distal-phalanx representations. Hence, we propose a 
well-ordered succession of the representations of the five digits 
as entities along the central sulcus, each entity comprised of the 
representations of the three phalanges. 

Within single digits, ordered phalanx representations (from 
pi to p3) along a similar direction (across subjects) could be 
shown for D5 and in a trend for D4. For the digits D3, D2, and 
Dl, the phalanx representations were neither generally ordered 
from pi to p3 nor were they aligned along a similar direc- 
tion across subjects, showing that the detailed arrangement of 
the phalanx representations within the intra-digit map varied 
between subjects. The lower degree of consistent intra-digit soma- 
totopy for these digits might be related to a higher degree of 
usage-dependent plasticity. 

This study underlines, based on high-resolution fMRI stud- 
ies of the entire digit area in human somatosensory cortex, the 
importance of quantitative across-digit and within-digit maps, 
to not only extract general organizational features, but also to 
adequately deal with the complexity of individual map structures. 
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